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1. Introduction 

In the continuous quest for Quark-Gluon Plasma (QGP), which began about two decades 
ago at the SPS at CERN, Geneva, and has continued in a spectacular way at the Relativistic 
Heavy Ion Collider (RHIC) at BNL, New York, USA, one now has a plethora of theoret- 
ical suggestions for signals of QGP as well as huge amounts of data. The Large Hadron 
Collider (LHC) at CERN will take us to even higher colliding energies of heavy ions than 
before. In the complex task of putting the pieces of the jigsaw puzzle together to establish 
eventually this new phase of strongly interacting matter from the experimental data, one 
needs as much of theoretical help as can be provided. In particular, theory needs to provide 
reliable estimates of various quantities, such as, the transition temperature, the equation of 
state (EoS) and the critical energy density needed to reach the QGP phase. Information 
on various properties of QGP, such as, the nature of its excitations and strength of their 
interactions would also be crucial in the experimental search. 

Quantum Chromo Dynamics (QCD) defined on a space-time lattice, lattice QCD, is the 
only successful and reliable tool to extract the desired non-perturbative physics from the 
underlying theory. This first-principles based and (essentially) parameter-free approach 
should be contrasted with the results from other approaches, primarily models. Thus not 
only does lattice QCD lead us to the phenomenon of quark confinement and spontaneous 
breaking of chiral symmetry (or the answer to why pion is so light) but it also provides 
us with a quantitative understanding of the spectrum of hadrons and their other proper- 
ties. Indeed, these results have been accorded an 'iconic' status in theoretical high energy 
physics recently [ 1 ] . One hopes that similar reliable information on the transition to QGP 
and properties of QGP will be provided by these techniques, and any experimental demon- 
stration of a failure of a prediction of lattice QCD, such as the transition to quark-gluon 
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plasma, will be tantamount to one of the best experimental evidence for physics beyond 
the standard model, discussed rather extensively at this workshop. 

Lattice formulation of QCD associates quark fields, ip(x), and the antiquark fields ip(x) 
with the site x = (x\,X2,Xs, x 4 ) of a 4-dimensional hypercubic lattice. The (inverse) lat- 
tice spacing a acts as the ultra-violet cut-off. Continuum limit of vanishing a corresponds 
to removal of the cut-off. As in the case of the continuum field theory, one obtains a lattice 
gauge theory by demanding invariance of the Lagrangian for free quark-antiquarks (e.g. 
obtained by a straightforward discretization of the usual Dirac Lagrangian) under any lo- 
cal phase rotation of these fields. This can be accomplished by introducing lattice gauge 
fields U£ = U fl (x) which are associated with a directed link from the site x to x + fia. 

Defining a partition function Z for these fields, akin to a complicated version of the 
familiar Ising model partition function, 



where Sg and Sf are gluonic and quark actions, M is the Dirac matrix in x, color, spin, 
flavour space for fermions of mass amj ea and a/j,f is the chemical potential (in lattice 
units). One can compute quantum expectation values of any physical observable O, which 
may contain fermion propagators of mass am valence , as averages with respect to the Z 
above. Thus, e.g., masses of physical particles are obtained from the exponential decays 
of appropriate correlation functions. Of relevance to heavy ion physics, is the investigation 
of eq. 1 at finite temperature, which is achieved by taking a N 3 x N t lattice, where N s is 
the number of lattice sites in a space direction and Nt in the time direction. This leads to 
volume V = N 3 .a 3 and temperature T = (aNt) -1 . Clearly, one needs to have N s » Nt. 

The Monte Carlo technique to evaluate the expectation value of any physical observ- 
able, consists of the following 3 steps : 1) Generate as large an ensemble of sets of links 
{Ug} for the whole lattice as possible, such that each set of {Ug} occurs with a proba- 
bility proportional to Ylf det Mf ■ exp[— Sg{{U£})], 2) Evaluate the observable for each 
configuration {Ug} and 3) Take its average over all the configurations in the set. Due to 
the enormity of the computational task to generate the set of {Ug} for full QCD, i.e., for 
a theory with all virtual quark loops included, one employs increasingly severe approx- 
imations with decreasing amount of computer time. These are (i) full QCD simulations 
on smaller lattices, (ii) partially quenched QCD simulations with am sea large and greater 
than am valence and (iii) quenched QCD simulations with am sea = oo (i.e. no dynamical 
quarks). The early lattice results and today's best results are obtained in the quenched ap- 
proximation. However, some aspects, such as the order of the phase transition, do depend 
strongly on the dynamical quark content, necessitating a judicious use of the quenched 
approximation. 

A transition temperature of T c ~ 170 MeV for two flavour QCD has been earlier es- 
timated [2,3]. The transition seems to be continuous. For 3 light flavours,the transition 
temperature seems to be lower by about 20 MeV [3], bracketing our world of two light and 
one heavy flavours, by that amount. Similarly, equation of state (EOS) has been predicted 
by lattice QCD. Other quantities, notably the Wroblewski Parameter A s , which measures 
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the strangeness enhancement in heavy ion physics, have also been predicted [4] by lattice 
QCD. These earlier results were obtained either in quenched approximation or on coarser 
N t =4 lattices for full theory. Thrust of the new results in the recent past has been on i) em- 
ploying larger lattices to achieve continuum limit and to lighter quarks, ii) more complex 
observables, such as speed of sound, transport coefficients, fluctuations and susceptibil- 
ities, J/i/)-dissolution/persistence, dileptons etc. and iii) theoretically more challenging 
T-/i phase diagram. An interesting comparison of predictions of conformally invariant 
QCD-like theories with lattice results has also been made. 

For reasons of both time and interest, I have chosen to limit this review to a few of the 
above mentioned selected topics. A quick overview of the basic lattice gauge theory can 
be found in [5] or many textbooks. In the next sections, I intend to discuss only the recent 
results, leaving out many technical details. A short summary is provided at the end. Let 
me emphasise here that one has witnessed recently a lot of activity in model building to 
explain the lattice QCD results. These include Quasiparticle models, Hadron Resonance 
Gas, Quarkonia from Lattice QQ potential, sQGP and coloured states. Many of them can, 
and eventually did, form interesting topics for working group discussions, leading even to 
a paper on the archive eventually. 

2. Properties of QGP 

In this section I shall review the recent progress made in the past few years in pinning 
down various properties of QGP using the lattice approach. These have been chosen for 
their direct connection with the experimental heavy ion physics. 

2.1 Equation of State and Speed of Sound 
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Figure 1. Recent results for EoS on Nt = 6 lattices from Refs. [6] (left) and [7] (right). 
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Equation of state plays a critical role in a variety of aspects of heavy ion collisions. 
Thus estimates of the energy density at T c are required to choose the appropriate heavy ion 
colliding energy. Recently results for the pressure and energy density for N t — 6 lattices 
at nonzero temperature QCD with 2+1 flavors of improved staggered quarks were reported 
[7]. Using improved gauge action and improved staggered quark action to eliminate the 
cut-off effects at 0(a 2 ) and with the heavy quark mass m s fixed at approximately the 
physical strange quark mass but with the two degenerate light quark masses m u d = 0.1m s 
or 0.2m s , it was found that the results were in agreement with those on Nt = 4 lattices, as 
seen on the right panel of Figure 1 . 

Using another improved action, similar results were obtained by another group [6] for 
quark masses corresponding to physical meson masses. These are displayed in the left 
panel of Figure 1 . Again, one sees small changes compared to the N t =4 case. In particular, 
the critical energy density, e(T c ) ~ 6T* still. One note of caution though is that the spatial 
volumes in physical units (such as inverse pion mass) are rather small and some changes 
may therefore be expected in true thermodynamic limit. 

Velocity of sound, C s , in QGP and dense hadronic media is a crucial input to many 
phenomenological studies based on the hydrodynamical approach. The presence of elliptic 
flow was established in this way and has been regarded as a key evidence in favour of 
the collective behaviour of the produced matter. Estimates of C s in the transition region 
have recently been obtained using a new approach [8] which relates it to the temperature 
derivative of the anomaly measure A /e, where A = e — 3P. Combining further with a 
new method to obtain e and P by an improved operator method, which leads to positive 
pressure on all lattices at all temperatures, C s as well as the specific heat C v were obtained 
[9] in the continuum limit. 
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Figure 2. C v /T 3 and 4e/T 4 (left) and C s 2 (right) as a function of T/T c for quenched 
QCD in continuum from Ref. [9]. 



Lattices with large temporal extent, Nt =8,10 and 12, and spatial sizes up to N s = 38 
were used in the quenched approximation to obtain the results shown in Figure 2. One 
sees in the left panel that C v ~ 4e for T > 2T C but its value is not close to the ideal gas 
limit. On the other hand, the speed of sound shown on the right, C 2 S , is close to the ideal 
gas limit by 2T C . Interestingly, it does not seem to show any structure near T c , whereas 
the specific heat does does hint at a peak at T c , It has been argued that fluctuations in px 
may be able to unravel the specific heat and its peak from heavy ion data. More theoretical 
studies are needed to sharpen this idea and to obtain a precise experimental measure. The 
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Figure 3. Ratio of shear viscosity and entropy density for quenched QCD [13]. 



above lattice results are encouraging for such studies. 

An exciting and curious agreement of the entropy density s computed for quenched 
continuum QCD has been reported [9] recently with that obtained using strong coupling 
prediction in the supersymmentric Yang-Mills theory in the temperature range of 2-3 T c . 
The latter predicts [10] s/s = f{g 2 N c ), where f(x) = § + §§C(3)£ _3/2 + • ■ • and 
s = ^n 2 N 2 T 3 is the ideal gas entropy density. In contrast, s = ^n 2 (N 2 — 1)T 3 for 
quenched QCD (i.e, without supersymmetry). Nevertheless, the prediction does seem to 
do well both in normalization and the shape in spite of the badly broken supersymmetry for 
the lattice results. Understanding this and extending it to even lower temperatures where 
weak coupling methods do fail would be very exciting indeed. 



2.2 QGP - (Almost) Perfect Liquid ? 

From the measured azimuthal angle distribution of various particles produced in a collision, 
one extracts a quantity called elliptic flow, V2- Phenomenological studies of ^(pt) show 
it to be consistent with QGP displaying ideal hydrodynamical behaviour [11]. Indeed, 
from the small deviations, one can derive [12] a bound on the shear viscosity, r\, of QGP: 
r s /7o = ^v/ s Tto < 0.1, where To is the formation time. This suggests a very small value 
for the dimensionless ratio, ij/s, whereas perturbation theory leads to a large value for it, 
giving rise to the description of QGP produced at RHIC as being an almost perfect and 
strongly coupled liquid. 

Kubo's linear response theory permits the determination of the transport coefficients, 
such as 77, in terms of equilibrium correlation functions. In particular, one needs the cor- 
relation function of the energy-momentum tensor. This is obtained on lattices in x-space 
and transformed to the momentum space, i.e, obtained at discrete Matsubara frequencies. 
In order to extract the transport coefficients, which are determined by the behaviour of the 
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retarded correlation functions at small frequencies u>, one continues these in the complex 
w-plane. Figure 3 displays the latest results [13] on the desired ratio ij/s on lattices as 
large as 24 3 x 8 but in QCD without dynamical quarks. Also shown are the estimates 
from perturbation theory and an analytic bound from supersymmetric QCD which is fairly 
close to the bound from hydrodynamics quoted above. The lattice results start from being 
close to these bounds in the vicinity of the transition and show a tendency of going towards 
the larger perturbative value as temperature increases. Calculations on larger lattices and 
with inclusion of dynamical quarks are needed to confirm these results and to make them 
more precise. Nevertheless their proximity to the bound from the heavy-ion data is very 
encouraging. 



2.3 Anomalous J/ij) Suppression. 
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Figure 4. Spectral functions of rj c , J/ip(left) and the \ states (right) in quenched 
QCD [16] at various temperatures indicated. 



Suppression of the most famous quarkonium, J/i/j, has been widely known as a signal 
for QGP production in heavy-ion collisions and a lot has been learnt [14] since the early 
data on this subject. Recently, the NA50 collaboration from CERN came out with their 
precise results. Moreover, J/?/>-suppression has been observed at the RHIC in BNL as well. 
To briefly summerise these latest results [15] from CERN and BNL and their comparison 
with model predictions, 

• Sulphur-Uranium and peripheral lead-lead (Pb-Pb) results for the ratio of J/ ip pro- 
duction cross section to that of Drell-Yan behave as expected from extrapolations of 
proton-nucleus data. Glauber model with an absorption cross section of 4.18 ± 0.35 
mb for the J/ip, determined from p-A data, decided these expectations. 

• Pb-Pb central data on the other hand show anomalous suppression with respect to 
the corresponding expectations. 

• The xjj'/DY ratio is not compatible for even S-U and peripheral Pb-Pb with similar 
expectations from proton-nucleus (with an absorption cross section of 7.6 ±1.12 
mb). 
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• ip' anomalous suppression sets in thus earlier than J/ip. 

• theoretical models that were successful in describing the above SPS, CERN data pre- 
dicted too much suppression at RHIC compared to what was observed by PHENIX 
atRHIC. 

The impressive sets of data from CERN and RHIC have attracted a lot of attention 
from lattice experts. The original Matsui-Satz idea of J/ip suppression was based on 
simple quarkonium potential models and an ansatz for the temperature dependence of the 
potential. It seems very important to check whether similar conclusions follows from the 
underlying theory, QCD, with as few assumptions as possible. Recognition of the so called 
MEM (Maximum Entropy Method) as a tool to extract the spectral functions of mesons 
from the temporal correlators computed on the lattice has permitted a fresh assessment of 
the original idea. Figure 4 exhibits the spectral functions [16] of the r} c , J/ip and the \ 
states in quenched QCD at various temperatures. 

As seen on the right panel, the peaks for Xc do not show any significance beyond the 
error bars by 1.1T C , i.e, they do seem to dissolve in agreement with the potential model 
estimate but the left panel shows persistence of J/ip and i] c up to 2.25 T c ; they seem to 
melt away only by 3T C . Similar results have been obtained by other groups [17], although 
differences persist on the precise melting temperature for J/ip. In addition to the need to 
iron them out, requiring larger lattices, uniform criteria and a comparison at the correlator 
levels, one has to also include dynamical quarks in these computations in quenched QCD. 
There is a lot more work to do but one can still ask whether these results should lead to 
changes in expectations of the suppression patterns as a function of temperature or the 
colliding energy since only a fifth to a third of the observed J/ip come from the states 
which do melt soon after the transition. Another interesting question is about quarkonia 
moving in the heat bath. One may expect them to see more energetic gluons, leading to 
more dissociation at the same temperature as the momentum increases. Preliminary results 
[18] show this to be indeed true for even J/ip and r\ c . However, the effect seems significant 
both below and above the transition, leading one to wonder whether it plays any role in the 
anomalous suppression seen in the heavy-ion collisions. 

3. QCD Phase Diagram 

Lattice QCD at nonzero baryon density may help us in understanding, or even deriving, 
an interesting physical phenomena, namely, color superconductivity, which may find ap- 
plications in the astrophysics of strange quark stars. From a theoretical viewpoint, it is, of 
course, crucial in completing the \ib~T phase diagram of QCD. Both the numerical and 
the analytical methods used at finite temperature are inadequate in this case due to the fact 
that the fermionic determinant det M (/j) is complex for /i ^ 0, commonly referred to 
as the sign (or the phase) problem. Various models, notably the Nambu-Jona-Lasinio or 
the random matrix model has played a big role in shaping our understanding of the QCD 
phase diagram. However, there have been some new exciting developments in the recent 
past for small fi. Most earlier attempts comprised of exploring first the zero temperature 
axis, where the problem is hardest. Recognizing this the latest strategy has been to work 
for small \i in the vicinity of the quark-hadron transition, and study its behaviour as /i is 
turned on. Various methods [19,20], such as the re-weighting method, a Taylor expansion 
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in ji, analytic continuation from imaginary fi, have lead to similar qualitative results. I will 
provide a flavour of these results by briefly mentioning here our results obtained by Taylor 
expansion, referring the reader to the original works in [19,20] for a detailed comparison. 

The Taylor expansion method has several advantages over the others. Prime amongst 
them is the ease of taking continuum and thermodynamic limit. As mentioned in Introduc- 
tion, one has to take these limits for the results to have any relevance to the real world of 
experiments. The re-weighting method, for example, has a factor that grows exponentially 
with lattice size and also has no systematic control over discretization errors. Analytic con- 
tinuation to real /i is also done term by term in a small /i-expansion. Employing the Taylor 
expansion, we [21] studied volume dependence at several T to i) bracket the critical region 
and then to ii) track its change as a function of volume. A strong volume dependence was 
found which changed the critical point obtained earlier by other methods substantially. The 
lattices we used were 4 xN%, with N s = 8, 10, 12, 16, 24, enabling us to vary the vol- 
ume V — N^a 3 at fixed temperature, i.e, fixed a. Our dynamical simulations for staggered 
fermions with two light dynamical (Nf = 2 of mass m/T c = 0.1) quarks were made using 
the well-known R-algorithm with a trajectory length scaled oc 7V S . From previous work 
[22], it is known that the transition temperature is T c /m p = 0.186 ± 0.006 and the choice 
of our quark mass corresponds to m^/m p — 0.31 ± 0.01. While still high compared to the 
real world, it is one of the smallest pion mass so far used; lowering pion mass further ne- 
cessitates even larger volumes than we were able to employ. Our simulations were made at 
T/T c = 0.75(2), 0.80(2), 0.85(1), 0.90(1), 0.95(1), 0.975(10), 1.00(1), 1.045(1), 1.15(1), 
1.25(2), 1.65(6) and 2.15(10). Typical statistics used was 50-100 in (max) autocorrelation 
units. 

Defining fif as the chemical potential for a flavour / = u,d,s and ^ = M« + + A*s 
and ^3 — ii u — iid as baryon and isospin chemical potentials, the respective density and 
susceptibility can be obtained from eq. 1 as : 

_ TdlnZ _ T d 2 lnZ 

ni ~ v~d^' Xij - (2) 

Setting fif = after taking the derivatives, n/ = but Xij are nontrivial. The diagonal 
X's are found [23] to be 

Xo = \pi{m u ) + ^0 2 {m u )] (3) 
X3 = \Oi{m u ) (4) 
Xs -J[Oi(m s ) + io 2 (m s )] (5) 

Here O t are trace of products of M~\ M' and M" and are estimated by a stochastic 
method: Tr A = £ti R\AR l /2N v , and (Tr A) 2 = 2£f >J=1 (Tr ^(Tr A) 3 /L{L - 
1), where Ri is a complex vector from a set of N v , subdivided in L independent sets. 
Further details can be found in [21,23]. 

Denoting higher order susceptibilities by Xn u ,n d , the pressure P has the expansion in /i: 

AP = Pjp, T) _ P(0,T) J_ /^\«- J_ (toy* 

T 4 ~ T 4 T 4 n u \\T) n d \\T) 

n u ,n d 

(6) 
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From this expansion, a series for baryonic susceptibility can be constructed. Its radius 
of convergence gives the nearest critical point. Successive estimates for the radius of con- 
vergence can be obtained from these terms using r n = , \ *+ 2 1 or p n = \\^w\] n ■ We 

y 1 x B 1 u Xb u 

used terms up to 8th order in /i, i.e., estimates from 2/4, 4/6 and 6/8 terms of the series eq. 
6. 

The ratio X11/X20 can be shown [21] to yield the ratio of widths of the measure in the 
imaginary and real directions at \i = 0. This argument for the measure of the imaginary 
part of the fermionic determinant can be generalized to nonzero /i with some care, by 
constructing similarly as above the coefficients for the off-diagonal susceptibility, xn{n)- 

Before going on further to the results on the critical point, let us take a detour of in- 
terest to the current heavy ion experiments. Quark number susceptibilities (QNS) which 
contribute the first nontrivial term in the above Taylor expansion also have their own in- 
dependent physical and theoretical relevance. They are crucial for some signatures of 
quark-gluon plasma such as fluctuations of charge and/or baryon number, and production 
of strangeness. Their additional theoretical importance is due to the check they provide on 
resummed perturbation expansions or any other scenario for the high T phase. 



3.1 The Wroblewski Parameter and Baryon-Strangeness Correlation 

Here we will touch upon two different aspects of the of physics hidden in the QNS. Ideally, 
one needs to obtain QNS in the continuum limit to extract any quantity of interest for heavy 
ion collisions. It has recently been shown that one can construct [24] robust ratios which 
have negligibly small theoretical and experimental systematic error. The primary amongst 
them is the Wroblewski parameter. Defined as the ratio of the strange particles and the 
non-strange particles produced in a collision, it has been studied widely as a measure of 
the strangeness production. Interestingly, most heavy ion collision data seems to yield a 
factor two higher value for it than other hadronic collisions. Using the continuum values for 
QNS, and under certain assumptions [4], one obtains X S (T C ) — 2Xs/(Xu + Xd) ~ 0.4 — 0.5 
which compares remarkably well with its latest RHIC value 0.47 ± 0.4, as shown in the left 
panel of Figure 5. While these results were obtained in quenched QCD in the continuum 
limit, also the full QCD results [24] shown in the right panel but on a small temporal lattice 
(N t = 4) are in very good agreement with these, as expected of a robust ratio. 

QNS can also be put to use in studying the correlation between quantum numbers K 
and L through the ratio C(kl)/l = ^(l'^-^ly^ = 2 ^ L - Again being ratios, these 
too are expected to be theoretically and experimentally robust. Such a strangeness-baryon 
number correlation, Cbs — ~ 3 ^r, was proposed [25] as a distinguishing test between 
the various models/pictures of quark-gluon plasma near T c . This, or the similarly defined 
strangeness-electric charge correlation, Cqs, is expected to be unity if quarks are the sole 
carriers of these quantum numbers (others being very heavy). On the other hand, the so- 
called sQGP model of Shuryak-Zahed [26] predicts these to be 0.66 and 1.2 respectively. 
Our lattice results [24] are exhibited in Figure 6. 

While one sees very different values for both Cbs an d Cqs below T c , their rapid ap- 
proach to unity clearly indicates the presence of quark-like degrees of freedom immediately 
above T c . In particular, the object which carries unit strangeness also has baryon number 
of -1/3 and a charge of 1/3, just like a strange antiquark would. We have found that a 
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Figure 5. Comparison of lattice results on the Wroblewski Parameter for quenched 
QCD with results from CERN and BNL (left). Corresponding full QCD results as a 
function of T/T c on the right. Taken from Refs. [4,24] respectively. 
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Figure 6. Css and Cqs as functions of T/T c . The quark masses used are 
Tn u d = 0.1T C and m s = T c . Taken from Ref. [24]. 



variation of the strange quark mass, m s /T c , between 0.1 and 1.0 does not alter either the 
value for T > T c , or its T-independence. A natural explanation of the T-behaviour arises 
if strange excitations with baryon number become lighter at T c . T-independence further 
suggests existence of a single such excitation. Such a picture is tantalizingly close to the 
canonical expectations of a transition from hadrons to quarks. 



3.2 The Critical Point 



Using both the definitions the radius of convergence above for the terms up to the 8th order 
in /i in pressure, (6th order in baryon number susceptibility), one obtains successively 
better estimates order by order. This can be done on each spatial volume, leading to a 
study of the volume dependence of the radius of convergence. Figure 7 shows the results 
for both p n and r n on our smallest (8 3 ) and largest (24 3 ) lattice. 

On the smaller lattice our results are consistent with the earlier estimate using re- 
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Figure 7. Radii of convergence p n (left) and r n (right) as a function of the order of 
the expansion at T — 
Taken from Ref. [21]. 



the expansion at T = 0.95T C on a 4 x 8 3 lattice (circles) and a 4 x 24 3 lattice (boxes) 



weighting method. We observe strong finite size effects around N a m v ~ 6 but our results 
on the largest lattice suggest a good stability with respect to an increase in the order of 
the expansion. Extrapolation in n leads to the following estimate for the critical point : 
fi E /T E = 1.1 ± 0.2 at T E = 0.95T C . Similar estimate has also been obtained for the 
re-weighting method recently although the quark mass used ( and hence the mass of the 
pion) was smaller in that case [27]. Another attempt to look [28] for the critical point using 
terms only up to 6th order did not find any critical point but they used much larger quark 
masses ( m^/irip ~ 0.7). 

1.1 i , , , , 1 
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H B /T 

Figure 8. QCD phase diagram and the freeze-out curve superimposed. The filled and 
open circles denote our [21] and earlier [19] estimate of the critical point for nearly the 
same quark mass. A scale of the CM energy per nucleon, yS, has been marked on the 
freeze-out curve. 
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4. Summary 

Lattice QCD predicts new states of strongly interacting matter and is able to shed light on 
the properties of the Quark-Gluon Plasma phase. One of the major developments of the 
recent past in lattice QCD is the firming up of the QCD phase diagram in \i-T plane on 
small N t . Different fermions different methods of simulations for nonzero /i, all lead to 
good agreement on the qualitative as well as the quantitative aspects. All estimates of T c , 
and (Te, are mutually consistent when compared for the right quark masses and in 
the thermodynamic limit. Our estimate for the critical point is hb/T ~ 1 — 2, as shown 
in Figure 8. Also shown in the figure is a freeze-out curve [29], converted to a value of 
T c appropriate to our computation. Such a curve results from the analysis of the heavy- 
ion data on particle yields. The required collision energy to reach the appropriate point 
are marked on the figure, which indicate the exciting possibility of discovering the critical 
point in a low energy RHIC run in the near future. 

Various physical quantities have been obtained in the continuum limit in the quenched 
approximation to QCD. These include the equation of state, the specific heat, the speed of 
sound in the neighbourhood of T c and the quark number susceptibilities. While the former 
are needed in hydrodynamical analysis of the particle spectra, and the resultant collective 
flow, the latter (QNS) are directly relevant to the physics of quark-gluon plasma signals 
at RHIC. The quenched lattice QCD estimate of Wroblewski parameter, A s , which is a 
measure of strangeness production in heavy ion collision experiments, has been known to 
be in excellent agreement with the RHIC and SPS results. 

Strong efforts are going on to extend these results to full QCD. Many features seem to 
change very little, although the major change in form of the order of the phase transition 
leads to quantitative differences near T c . Interestingly the robust Wroblewski parameter 
seems also not to change quantitatively near T c , although one still needs more precise 
results and that too for realistic pion and kaon masses. Our results on baryon number- 
strangeness and electric charge-strangeness correlations suggest i) a rapid change in going 
through the transition in full (Nf = 2) QCD, and ii) the quark-gluon plasma to have quark- 
like excitations even close to T c . 

The heavy ion data from CERN and BNL has provided us a lot of surprises, and will 
continue to do so in future. Lattice QCD has played an important part in understanding 
some of these, and clearly a lot more work is ahead both in form of better precision as well 
as new ideas for experimentally measurable lattice predictions. 
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